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Context. White dwarf evolution is essentially a gravothermal cooling process, which, for cool white dwarfs, sensitively depends on 

■ the treatment of the outer boundary conditions. 
Aims. We provide detailed outer boundary conditions appropriate for computing the evolution of cool white dwarfs employing detailed 
non-gray model atmospheres for pure hydrogen composition. We also explore the impact on the white dwarf cooling times of different 

■ assumptions for energy transfer in the atmosphere of cool white dwarfs. 
O [ Methods. Detailed non-gray model atmospheres are computed taken into account non-ideal effects in the gas equation of state and 

chemical equilibrium, collision-induced absorption from molecules, and the Lyman a quasi-molecular opacity. 
Results. We explore the impact of outer boundary conditions provided by updated model atmospheres on the cooling times of 0.60 and 
0.90 Mq white dwarf sequences. Our results show that the use of detailed outer boundary conditions becomes relevant for effective 
temperatures lower than 5800 and 6100 K for sequences with 0.60 and 0.90 Mq, respectively. Detailed model atmospheres predict 
ages that are up to « 10% shorter at log(L/LQ) = -4 when compared with the ages derived using Eddington-like approximations 
at Tr oss = 2/3. We also analyze the effects of various assumptions and physical processes of relevance in the calculation of outer 
boundary conditions. In particular, we find that the Lyce red wing absorption does not affect substantially the evolution of white 
dwarfs. 

1 Conclusions. White dwarf cooling timescales are sensitive to the surface boundary conditions for r e ff ~ 6000 K. Interestingly enough, 

ly-^ , non-gray effects have little consequences on these cooling times at observable luminosities. In fact, collision-induced absorption 

. processes, which significantly affect the spectra and colors of old white dwarfs with hydrogen-rich atmospheres, have not noticeable 

■ effects in their cooling rates, except throughout the Rosseland mean opacity. 

■ Key words, stars: evolution — stars: interiors — stars: white dwarfs 

o: 

I \ 1. Introduction (Mestel 1952). Once convection reaches the outer edge of the 

J> . degenerate core in low-luminosity white dwarfs - the so-called 

• r-j ■ An accurate assessment of the rate at which white dwarfs cool convective coupling - the cooling becomes strongly tied to the 

P>^ ; down is a fundamental issue, because these stars can be used treatment of the outer boundary conditions (Bohm & Grenfell 

H . as independent and accurate age indicators. As a matter of fact, 1973). Consequently, an accurate assessment of the cooling rate 

white dwarfs are the most common end-point of stellar evolu- at low luminosities requires the use of detailed model atmo- 

tion - see, for instance, Althaus et al. (2010a) for a recent re- spheres (Hansen 1998; Salaris et al. 2000; Serenelli et al. 2001). 
view - and as such are valuable in constraining several properties 

of a wide variety of stellar populations, including globular and The treatment of the energy transfer in white dwarf atmo- 
open clusters (Von Hippel & Gilmore 2000; Hansen et al. 2007; spheres is a difficult task which involves the solution of the equa- 
Winget et al. 2009; Garcfa-Berro et al. 2010). Additionally, they tions of radiative transfer coupled to convection, in a highly non- 
can be used to place constraints on elementary particles such as ideal gas regime where several molecular and quasi-molecular 
axions (Isern et al. 1992; Corsico et al. 2001; Isern et al. 2008), processes have to be considered. This implies a high degree of 
and neutrinos (Winget et al. 2004) or on alternative theories of sophistication of the calculations, especially at very low lumi- 
gravitation (Garcfa-Berro et al. 1995; Garcfa-Berro et al. 201 1). nosities. The importance of using detailed boundary conditions 
These and other potential applications require a detailed and pre- was first addressed by Hansen (1998, 1999). Over the years, de- 
cise knowledge of the main physical processes that control their tailed model atmospheres that include a complete treatment of 
evolution. A key ingredient is the energy transfer in the atmo- energy absorption processes such as collision-induced opacity 
spheric and sub-atmospheric layers which control their cooling (Bergeron et al. 1991; Saumon and Jacobson 1999; Rohrmann 
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Fig. 1. Test of the diffusion approximation in the plane (T e g, 
logTR OSS ) based on differences between the ratio K v /J v and its 
asymptotic value 1/3 at large optical depth, for a wavelength 
A = 2000 A. These results correspond to detailed non-gray mod- 
els for H atmospheres at log g = 8. 



2001) and the Lyman a wing absorption (Kowalski & Saumon 
2006; Rohrmann et al. 201 1), have been developed. 

Here we provide detailed boundary conditions which allow 
to consistently compute the evolution of cool white dwarfs with 
pure hydrogen atmospheres. These boundary conditions are pro- 
vided in the form of tables for a wide range of surface gravities 
and effective temperatures. In the following sections we describe 
the model atmospheres and the evolutionary code (§|2j, and ex- 
plore the impact on the cooling times of different physical pro- 
cesses (§(3). Conclusions are given in Sec. [4] 

2. Numerical tools 

In our calculations, the outer boundary conditions are obtained 
using the pure-hydrogen LTE model atmospheres described at 
length in Rohrmann et al. (2001, 2002, 2011). Specifically, we 
compute pressure, temperature, and outer mass fraction at an 
Rosseland mean optical depth tr oss = 25.1189 (logTR OSS = 1.4) 
for 40000 K < T eS < 2000 K and 6.5 < logg < 9.5. Model 
atmospheres were computed in the range -6 < log tr oss < 2 (in 
steps of 0. 1 dex) assuming hydrostatic and radiative-convective 
equilibrium. Convective transport is treated within the usual 
mixing-length (ML2) approximation, in which the ratio of the 
mixing-length to the pressure scale height is a — I. The micro- 
physics comprises non-ideal effects in the gas equation of state 
and chemical equilibrium based on the occupation probability 
formalism as describe in Rohrmann et al. (2002). The chemi- 
cal composition of the atmosphere includes H, H2, H + , H , H^, 
H+, He, He , He\ He 2+ , He+, HeH + , and e . The level occu- 
pation probabilities are self-consistently incorporated in the cal- 
culation of the line and continuum opacities. Collision-induced 
absorptions due to H2-H2 (Borysow, Jorgensen & Fu 2001), 
H 2 -H (Gustafsson & Frommhold 2003), H-H (Doyle 1968), 
H2-He (Jorgensen et al. 2000) and H-He pairs (Gustafsson & 
Frommhold 2001) are also taken into account. Model atmo- 
spheres explicitly include the Lyman a quasi-molecular opac- 
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Fig. 2. Same as Fig. [Q but for A = 5000A. 



ity as a result of eight allowed electric dipole transitions aris- 
ing from H-H and H-H2 collisions (Rohrmann et al. 201 1). This 
opacity reduces the predicted flux at wavelength A < 4000A for 
stars cooler than T e ff « 6000 K. 



3. Results 

The evolutionary calculations reported here have been done us- 
ing the LPCODE stellar evolutionary code (Althaus et al. 2003, 
2005, 2012). This code has recently been successfully used to 
perform very accurate evolutionary calculations - see Garcfa- 
Berro et al. (2010), Althaus et al. (2010b), Renedo et al. (2010), 
and references therein. Of relevance for this work, outer bound- 
ary condition to the stellar structure and evolution equations are 
specified by performing three envelope integrations from start- 
ing values - as given by the adopted model atmosphere - inward 
to a fitting outer mass fraction, as described in Kippenhahn et al. 
(1967). Energy sources resulting from crystallization - the re- 
lease of latent heat and of energy resulting from carbon-oxygen 
phase separation - are taken into account using the phase di- 
agram of Horowitz et al. (2010), see Althaus et al. (2012) for 
details. The equation of state is that of Segretain et al. (1994) 
for the high-density regime - which accounts for all the impor- 
tant contributions for both the liquid and solid phases (Althaus 
et al. 2007) - complemented with an updated version of the 
equation of state of Magni & Mazzitelli (1979) for the low- 
density regime. Radiative opacities are those of OPAL (Iglesias 
& Rogers 1996), including carbon- and oxygen-rich compo- 
sitions, complemented with the low-temperature opacities of 
Ferguson et al. (2005), linearly extrapolated to high densities 
when needed. Conductive opacities are taken from Cassisi et al. 
(2007). We also take into account the effects of element diffusion 
due to gravitational settling, chemical and thermal diffusion of 
'H, 3 He, 4 He, 12 C, 13 C, 14 N and 16 0, see Althaus et al. (2003) for 
details. In particular, the metal mass fraction Z in the envelope 
of our models is specified by scaling it to the local abundance of 
the CNO elements at each layer. To account for this, we consider 
radiative opacities tables from OPAL for arbitrary metallicities. 
Convection is treated within the ML2 version (a = 1) of the 
mixing-length theory. 
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Fig. 3. Fraction of convective energy in the plane (log tr oss , T e g) 
for model atmospheres with \ogg = 8. 



We compute the evolution of sequences of white dwarfs 
of 0.6 and 0.9 Mq. Initial configurations are the result of the 
complete evolution of 1.75 and 5.0 Mq progenitors, respec- 
tively, with metallicity Z = 0.01 - see Renedo et al. (2010) 
for details. Progenitor stars were evolved from the zero age 
main sequence, through the thermally-pulsing and mass-loss 
phases on the asymptotic giant branch, to the cooling phase. 
Time-dependent overshoot mixing beyond the formal convective 
boundary during core hydrogen and helium burning stages was 
taken into account, see Althaus et al. (2005) for details. Mass 
loss during the RGB and AGB phases was considered following 
the prescription of Schroder & Cuntz (2005) and Vassiliadis & 
Wood ( 1 993), respectively. The outer chemical profiles of our se- 
quences are the result of element diffusion processes that lead to 
the formation of pure hydrogen envelopes with zero metallicity. 
The total mass of hydrogen left after hydrogen burning is com- 
pletely exhausted amounts to 7 x 10~ 5 and 7.6 x 10~ 6 Mq for the 
0.6 and 0.9 Mq sequences, respectively. For each stellar mass, 
we computed the cooling phase down to very low luminosities, 
when most of the white dwarf has already crystallized. 

Standard outer boundary conditions are usually based on the 
Eddington gray approximation, which assumes the diffusion ap- 
proximation for radiative transfer and neglects convection at low 
optical depths lower than tr oss = 2/3. A simple test for the dif- 
fusion approximation is to compare the ratio between the second 
(K v ) and zero-order (J v ) moments of the radiation field, with its 
asymptotic value (1/3) at large depth (tr oss — > oo). Figs. Q] and [2] 
show such comparison at optical depths 10~ 4 2 < tr oss < 10 2 
for detailed, pure hydrogen atmospheres with \ogg - 8 and 
2000 K < r eff < 19000 K, and for A = 2000 and 5000A, 
respectively. In general, the difference between K v /J v and 1/3 
increases outward, first with negative values and then with posi- 
tive ones, forming a complex pattern in the plane (T e g, log tr oss ). 
Departures from the diffusion limit at Tr oss = 2/3 (logTR OSS = 
-0.176) become larger than l%-5% depending on the wave- 
length and effective temperature. Large deviations occur mainly 
for the visible (Fig. [2j and infrared wavelengths at low T e g. This 
depth, therefore, is not optimal to establish the outer boundary 
conditions of cool white dwarfs. 



Convective transport represents another serious limitation of 
the standard method to evaluate boundary conditions. Fig. [3] 
shows the fraction of the energy flux carried out by convection 
as a function of T e g in model atmospheres with \ogg — 8. A 
superficial convection zone starts at T e ff « 18000 K associated 
with the recombination of hydrogen. Below T e g a 12200 K, 
the efficieny of convection increases rapidly with depth and the 
convection zone extends down to the bottom of the atmosphere 
(tross = 100). When the star cools below r eff a 7000 K, the 
top of the convection zone slowly extends to very low optical 
depths as a result of H2 formation. Convection efficiency de- 
clines for models cooler than r e g » 4000 K. It is clear that for 
3000 K ;$ T e ff S5 12500 K, convection may carry more than 
30^0% of the total flux at r Ross < 2/3. 

The importance of the gray approximation throughout the 
atmosphere may be tested comparing two mean opacities with 
different spectral weighting, for instance the Planck (/fpianck) and 
the Rosseland (a-r oss ) means (Mihalas 1978). Note that all mean 
opacities are equal in a gray atmosphere. Fig. [4] displays Kpi an ck 
and /fR 0SS for a hydrogen gas as a function of the temperature, for 
several densities. At high (log T ^ 4) and low (log T £ 3.5) tem- 
peratures, the opacity is mainly due to atomic and molecular hy- 
drogen, respectively, both yielding large discrepancies between 
Kpianck and /cross, and therefore strong deviations from the gray 
approximation. At intermediate temperatures (3.5 5= logT ^ 4), 
the differences between the two mean opacities are smaller be- 
cause the H absorption, which is the dominant opacity source, 
has a nearly flat behavior in part of the spectrum. Fig. [4] also 
shows the run of Kpi an ck (dashed lines) and /fR 0SS (dash-dotted 
lines) for H atmospheres with logg = 8 and r e ff = 3000 K, 
8000 K, and 30000 K (from left to right). As expected, the 
model with T e g = 8000 K shows moderate differences be- 
tween both opacities (ss 0.5 dex). Discrepancies in the hot model 
(T eS = 30000 K) are somewhat larger (0.6 - 1.2 dex), but the 
largest ones (1-3 dex) occur in the cooler model (T e g = 3000 
K), which strongly deviates from the gray approximation. 

The effects of quasi-molecular processes on the Rosseland 
mean opacity are illustrated in Fig.[5]for logp = -8 and — 1, The 
dashed and dashed-dotted lines correspond respectively to cal- 
culations in which the contributions of the Lya red wing and the 
H2-H2 CIA opacities have been removed. While the Lya opacity 
affects the energy distribution emitted by white dwarfs cooler 
than r e ff ~ 5000 K (Kowalski & Saumon 2006; Rohrmann et 
al. 2011), it has a moderate effect on the Rosseland opacity at 
low densities and relatively high temperatures (see Fig.|5]l, when 
the convective coupling with stellar core has not yet occurred. In 
contrast, abrupt changes in kr oss occur at low temperatures due 
to the H2-H2 CIA processes, which increase the mean opacity 
by several orders of magnitude and have important effects on the 
cooling rates of old white dwarfs (Hansen 1999). 

In view of the previous remarks, it seems important to treat 
the boundary conditions as accurately as possible. Detailed non- 
gray model atmospheres have been computed to provide surface 
boundary conditions for white dwarfs with hydrogen envelopes. 
The values of the pressure, temperature and outer fraction of stel- 
lar mass were obtained at tr oss = 25.1189 (logTR OSS = 1.4) for 
effective temperatures ranging from 2000 to 40000 K, in steps 
of 100 K, and \ogg from 6.5 to 9.5, in steps of 0.1 dex. Note that 
at tr ss > 25, the diffusion approximation is guaranteed within 
0.01% or better for most of the spectrum (see Figs.Q]and[2]i. 

To explore the influence of the boundary conditions on the 
cooling, we display in Figs. [6] and [7] the relationship between the 
surface luminosity and age for the 0.6 and 0.9 Mq sequences, 
respectively, that result from non-gray model atmospheres (solid 
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Fig. 4. Planck (solid lines) and Rosseland (dotted lines) mean 
opacities as a function of the temperature at mass densities in- 
creasing from logp = -10 to logp = -1 in steps of one dex 
(from bottom to top). Thick dashed and dashed-dotted curves 
are the Planck and Rosseland mean opacities for white dwarf 
atmospheres with logg = 8 and r eff = 3000 K (left), 8000 K 
(middle) and 30000 K (right). 



line), gray model atmosphere (dotted-dashed line) and gray 
model atmospheres in which the convection is neglected (dotted 
line). These figures also show results obtained with boundary 
conditions at tr oss =2/3 based on the Eddington gray approxi- 
mation (dashed line). Clearly, the use of detailed outer boundary 
conditions becomes relevant for cooling once evolution has pro- 
ceeded to luminosities lower than log(L/Lo) ~ ~3.8 (-4.0) in 
the case of 0.60 Mq (0.90 Mq) models. These luminosities cor- 
respond to effective temperatures lower than 5800 K (6100 K). 
Such values indicate the onset of the convective coupling be- 
tween the outer envelope and the isothermal degenerate core 
(Tassoul et al. 1990; D' Antona & Mazzitelli 1990; Prada Moroni 
& Straniero 2007). For larger effective temperatures, evolution is 
almost insensitive to a detailed treatment of the outer boundary 
conditions. 

Figs. [6] and [7] show that the use of detailed model atmo- 
spheres directly translates into cooling times different from those 
predicted by the standard Eddington approximation. Note that, 
for the 0.6 Mq sequence, in the range -3.8 <; log(L/Lg) ~ -4.3, 
the inclusion of proper outer boundary conditions decreases the 
cooling ages by up to 0.7 Gyr, while this trend in the cooling 

times is reversed at luminosities below log(L/L0) 4.3, where 

the use of detailed model atmospheres results in longer cooling 
times. The behavior of the cooling times is qualitatively similar 
for the more massive sequence. Again, the use of any Eddington- 
like approximation that involves the diffusion assumption for ra- 
diative transfer and neglects convection at low optical depths, 
incorrectly predicts the evolution of cool white dwarfs. 

Surprisingly, and contrary to what was expected from Fig. [4] 
Figs. [6] and [7] show no appreciable changes in the cooling times 
when the gray approximation is assumed in the model atmo- 
spheres. In these calculations the monochromatic opacity co- 
efficient was forced to take the value of the Rosseland mean 
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Fig. 5. Effects of the collision-induced opacities (H2-H2 CIA and 
Lya red wing) on the Rosseland mean opacity for two densities, 
logp = -1 and logp = -8. 



opacity. In fact, Fig.[8]shows that at the onset of convective cou- 
pling (T e ff ~ 6000 K) the temperature-pressure stratifications 
of gray and non-gray models are practically identical. Below 
T e g x 5000 K, H2-H2 collision-induced opacity reduces the 
surface temperature and increases the temperature in deep at- 
mospheric layers respect to the gray model. These non-gray ef- 
fects increase towards r e ff ~ 3000 K (Fig. [8]), which corre- 
sponds to about log(L/L0) ~ -4.95 for the 0.6 Mq sequence 
(see Fig. [6]), but nevertheless have little consequences on the 
cooling times at observable luminosities. Indeed, they are promi- 
nent at very low T e s, and are expected to influence the cooling 
times once evolution has proceeded to luminosities lower than 
log(L/L0) ~ -4.95. We have also computed cooling sequences 
where the collision-induced broadening of Lya was omitted in 
model atmospheres, and we found that this opacity source does 
not affect the evolution substantially. These results show that 
processes which markedly alter the distribution of spectral en- 
ergy radiated by the star may have no effect on its cooling time. 

On the contrary, neglecting convective energy transfer in 
model atmospheres strongly alters the cooling times. As illus- 
trated in Fig. [8] neglecting convection results in much larger 
temperatures at the base of the atmosphere (dotted lines in the 
figure), thus producing a markedly shallower outer convection 
zone, and eventually resulting in age differences of up to 2 Gyr. 
In particular, for a 0.6 Mq model at log(L/LQ) = -4, the use of 
outer boundary conditions derived from model atmospheres that 
neglect convection down to an optical depth of tr oss as 25 leads 

to an outer convective zone with mass log(M conv /MwD) 18, 

whereas the value resulting when convection is considered is 
-13.5. 

The impact on the cooling times of the different bound- 
ary conditions can be better appreciated by inspecting Fig. [9] 
which illustrates the age differences with respect to the case in 
which the Eddington gray approximation is used for sequences 
that consider non-gray model atmospheres (solid blue line), gray 
model atmospheres (dotted-dashed green line), and gray model 
atmospheres in which convection is neglected (dotted line). The 



4 



R. D. Rohrmann et al.: Boundary conditions for evolving cool white dwarfs (RN) 



-3.6 



-4.8 



1 

V 


I 1 
m wd = 0-60 M sun 


- 


•"V N 
• x v 




- 








Eddington (t=2/3) 

gray 

gray (no convection) 

non-gray 


\ \. .. 
\ \ 

\ \ 
\ V 
s V 
\ \ - 
\ \ 
\ \ 






\ \ 
\ \ 
\ \ 
\ \ 
\ \ 

, i , W 



4 6 8 10 12 

Age (Gyr) 



14 



Fig. 6. Surface luminosity versus age for the 0.6 Mq cooling 
sequences resulting from the use of outer boundary conditions 
as given by non-gray model atmospheres (solid blue line), gray 
model atmospheres (dotted-dashed green line), and gray model 
atmospheres in which convection is neglected (dotted line). 
Results based on the Eddington gray approximation are shown 
with a dashed red line. 
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ferent assumptions as indicated on the plot. The top and bottom 
layers of each model are located at tr oss = 10~ 6 and tr oss = 100, 
respectively. The thick lines show the layers located between 
tross ~ 2/3 and t Ross « 25. 
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Fig. 7. Same as in Fig.|6]for the 0.9 Mq sequences. 



upper and bottom panel correspond to the 0.6 and 0.9 Mq se- 
quences, respectively. Note that for the 0.6 Mq sequence age 
differences are negative in the range -3.8 <: log(L/LQ) ^ -4.3, 
and the use of detailed model atmospheres predicts ages that are 
up to * 10% shorter at log(L/LQ) = -4.05 when compared with 
the ages derived using the Eddington approximation. The differ- 
ences are somewhat smaller for the 0.9 Mq sequence, reaching 
up to « 7% at log(L/LQ) = -4.2. Besides, the gray assumption 
translates into age differences less than ~ 1% in both sequences. 



Since the input physics adopted in the codes used to compute 
the stellar interior and the atmosphere is not exactly the same, 
we have examined the impact of matching atmosphere and inte- 
rior models at different optical depths. We found that the cool- 
ing times differ at most 0.7% (0.9%) for boundary conditions at 
tross ~ 50 (100) respect to those obtained at tr oss w 25. Thus, 
differences in the constitutive physics in the codes appear to have 
small consequences in the derivation of the boundary conditions. 
We have also tested the effects of changing the value of the mix- 
ing length parameter (a) in the convection theory. In particular, if 
the efficiency of convection is increased to a — 2 (ML3 version 
of the mixing-lenght theory), the relative age differences respect 
to the use of a = 1 (ML2 version) becomes smaller than 0.8% 
(0.4%) for the 0.6 Mq (0.9 Mq) model. 

4. Conclusions 

The purpose of this work has been to provide detailed outer 
boundary conditions which allow to compute white dwarf evo- 
lution in a consistent way with the predictions of detailed model 
atmospheres. Data are provided in the form of tables for a wide 
range of surface gravities and effective temperatures, which are 
appropriate for computing the evolution of cool white dwarfs 
with pure hydrogen atmospheres. The full set of data is avail- 



able at http://www.icate-conicet.gob.ar/rohrmann/tables.html or 
upon request to the authors at their e-mail addresses. 

White dwarf cooling timescales are sensitive to the sur- 
face boundary conditions for effective temperatures lower than 
T e g w 6000 K. Different outer boundary conditions may result 
in substantial differences in the cooling times for cooler white 
dwarfs. However, non-gray effects do not become important in 
the cooling rates. On the other hand, depending on the stellar lu- 
minosity, the use of detailed model atmospheres like the ones 
presented here results in age differences of about 10% when 
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Fig. 9. Age differences between sequences that consider non- 
gray model atmospheres (solid blue line), gray model atmo- 
spheres (dotted-dashed green line), and gray model atmospheres 
in which convection is neglected (dotted line) with respect to the 
case in which the Eddington gray approximation is considered. 
The upper and bottom panel correspond, respectively, to the 0.6 
and 0.9 Mq white dwarf sequences. 

compared with the ages computed using the Eddington approxi- 
mation, which assumes the diffusion approximation for radiative 
transfer and neglects convection at low optical depths lower than 
tross =2/3. These differences are on the order of the current un- 
certainties in the white dwarf cooling times at low luminosities 
that result from uncertainties in the treatment of progenitor evo- 
lution, particularly during the core helium burning phase (Prada 
Moroni & Straniero 2002; Salaris et al. 2010). Consequently, 
accurate outer boundary conditions provided by detailed model 
atmospheres have to be considered in evolutionary studies aimed 
at using these stars as accurate cosmic clocks. 
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